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Abstract 

We develop a mode coupling theory(MCT) to study the nonequilibrium glass transition behav¬ 
ior of a mono-disperse mixture of active-passive hard-sphere particles. The MCT equations clearly 
demonstrate that the glass transition is shifted to higher values of total volume fraction when 
doping a passive system with active particles. Interestingly, we find that the glass transition point 
may show a non-monotonic dependence on the effective diffusivity of the active component, indi¬ 
cating a nontrivial type of activity induced reentrance behavior. Analysis based on the nonergodic 
parameters suggest that the glassy state at small activity is due to the caging effect, while that at 
high activity could result from activity induced dynamic clustering. 
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I. INTRODUCTION 


The collective behavior of systems containing active particles have gained extensive atten¬ 
tion in recent years due to its great importance both from a fundamental physics perspective 
and for understanding many biological systems [l, 2|. A wealth of new nonequilibrium phe¬ 
nomena have been reported, such as active swarming, large scale vortex formation js, 4|, 
phase separation js-Q], etc, both experimentally and theoretically. Recently, a new trend in 
this field has been the dynamics of dense assemblies of self-propelled particles around glass 
transition. Experiments on crowded systems of active colloids and active cells show dy¬ 
namic features such as lamming and dynamic arrest that are very similar to those observed 
in glassy materials [lO, 11|. Computer simulations demonstrated that nonequilibrium glass 
transition or dynamic arrest behavior does occur in a dense suspension of self-propelled hard 
[here the critical density for glass transition shifts to larger value with increasing 
13| . Starting from a generalized Langevin equation with colored non-thermal 
noise, L.Breather and J. Curtain theoretically predicted that dynamic arrest can occur in 
systems that are far from equilibrium, showing that non-equilibrium glass transition moves 
to lower temperature with increasing activity and to higher temperature with increasing 
dissipation in spin glasses 1^. Mode coupling theory (MCT) were also proposed recently to 
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study the glassy dynamics of driven granular fluids 


15| and active colloidal suspensions jl6|. 


While most of the studies so far have only considered single-component active particles, 
very recently, mixture ^sterns of active-passive particles began to draw new attentions. 
Interesting experiments p, ll7| reported that the addition of active particles in a system 
can dramatically alter its phase behavior. Molecular dynamics simulation showed that 
introduction of activity to a passive system may not only hamper phase separation, but can 
enhance it as well, based on the coordination among the active particles 18|. It was also 
demonstrated th^ activity can induce phase separation and direct self-assembly in active- 

ations showed that one may crystallize 


passive mixtures}^. In particular. Brownian simu! 
hard-sphere glass by doping with active particles 


191], which proposes a very interesting 


question about the glass transition behavior of active-passive mixtures. However, an unihed 
microscopic theoretical framework to describe this important issue is still lacking. 

In the present paper, we develop a general MCT framework to study the glass transition 
behavior of active-passive mixtures. Our starting point is the Smoluchowski equation for N- 
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particle probability density function, wherein the particle activity is realized via an effective 
diffusivity that is larger than that of passive particles. Such a treatment allows one to apply 
the Mori-Zwanzig projection operator formalism such that a set of closed equations regarding 
the time evolutions of the density correlators can be obtained. In particular, we apply this 
approach to a binary mixture of mono-disperse active-passive hard sphere particles, with 
particular attention paid on how the particle activity and the number fraction of active 
particles would influence the glass transition behaviors. While doping with active particles 
can shift the transition to higher total volume fraction as expected, we find that increasing 
particle activity may lead to an nontrivial type of reentrance glass transition behavior. 


II. THEORY 


For generality, we consider an m-component mixture of N spherical colloidal particles, 
being active or not, dispersed in a simple fluid with temperature T and volume V. For passive 
particles, the over-damped dynamics can be described by Langevin equations involving 
the time evolution of the position vector rf for an /i-type particle labeled /. For active 
particles, one generally needs to consider a further orientation variable to account for 
the tumbling or rotational diffusion 20|, l2l| . Very recently, it was demonstrated that the 
dynamics of an active particle undergoing self-propulsion and rotational diffusion can be well- 
approximated by a random-walk in the long time limit with an effective diffusivity 
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These observations facilitate us to propose a minimal model for active-passive mixtures, 
which starts from the Smoluchowski equation for the probability density P (r'^, t) of the 
particle configuration 


dP{Y^,t) ldt = hP{T^,t) ( 1 ) 

where is the Smoluchowski operator 

m 

n = ( 2 ) 

fj.=l j=l 


Here denotes the bare diffusivity for /i-type particles with total number and (3 = 
l/ksT with ks the Boltzmann constant. = —V^f/ (r-^) is the direct force acting on 
the particle j G /i due to the total potential energy U (r^) of the colloidal particles, and 
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Vj is the gradient operator with respect to r^. A key factor here is the parameter 
characterizing the particle activity of /x-species: > 1 for active particles while A^ = 1 for 

passive ones. For an equilibrium distribution oc , '^jPe = + {P^j) Pe such 

that the random force balances the potential force for passive particles, but not for active 
ones. 

To probe the collective dynamics of the system, one generally considers the density cor¬ 
relator 


(q, t) = 


PH 


Po]P- 




( 3 ) 


where is the Fourier transform of the density p^(r, t) = 

—r^), and (■) denotes an equilibrium average, is the adjoint or back¬ 


ward Smoluchowski operator given by 
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Na 


fit = ^ ^ (A'^v; + /3F;) V 


( 4 ) 




i=i 


where for arbitrary functions / and g of we have J dr^ f ^fif?^ = f dr^ 


of adjoint Smoluc 
coupling methods 


lowski operator facilitates us to apply Mori-Zwanzig approach and mode- 


24i | to obtain an approximate dynamic equation for the density correlators 


which reads in the matrix form as (See the Supplemental Information) 


^$(g,t)+W(g)-S-'(g)-$(g,t) 

A- [ duM. (q, t — u) ■ iq) ■ (g, m) = 0 (5) 

Jo 

Herein, (g) is the inverse of the partial structure factor matrix S (g) = $(g,t = 0), 
W (g) is a frequency matrix given by 


W,, (g) = g^fio^ [(A^ - 1) S,, (g) + 5 ,,]. 


( 6 ) 


M (g, t) denotes the matrix of so-called irreducible memory function 24l-l26l| with elements 
given by 


(g, t) 


^ lit 

-2I: E (Q; k) (H’ k) 

k (57(5'7' = 1 

x$ 55 / (fc,t) (|q - k| ,t) 


( 7 ) 
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with vortex functions 


Vf.-S'r (q, k) 


-^|(q.k)i„C„nk) 

V 

+q ■ (q - k) (|q - k|)] 


( 8 ) 


where p = N/V is the total number density and {k) = {k) in which c^y {k) 

are direct correlation functions. C^y {k) is related to the static structure factor via pC^y = 

V - . 

The MCT equations ([5]) to ([8]) for general active-passive mixtures constitute the hrst part 
of central results of the present paper. In general, one may employ them to study the glass 
transition or dynamic arrest behavior of any multi-component systems. In the present work, 
we will mainly focus on a two-component hard-sphere system, one is active labeled by ’a’ 
and the other is passive labeled by ’p’ . For simplicity, we consider that the diameters d 
of both types of particles are the same. The total volume fraction is given by = p7rd^/6, 
with d set to be 1. The number fraction and the activity parameter of the active component 
are given by Xa and Xa respectively, which are chosen as the main control parameters in the 
present study. Eqs.© are numerically solved with Percus-Yevick static structure factors 
S^iu {k) = a or p) as input and setting Dq = = 1. For a pure passive system with 

Xa = 0, the equations predict a glass transition (GT) at volume fraction = 0.515. 


III. RESULTS AND DISCUSSION 

First of all, we study how the particle activity influences the GT. To this end, we £x 
Xa = 0.5 and set p = 0.524 which is above the GT point for pure passive system. 
In Fig{T](a), we show the time evolution of the dynamic scattering functions ^aa {k,t) and 
<hpp {k,t) evaluated at /c = km, where km — 7.23 is the location of the hrst peak in Sk, with 
varying particle activity A“ in a relatively small range. When A“ is just slightly than one, 
e.g., A“ = 1.02, both correlators do not decay in the long time and the system remain in 
the glassy state. The nonegodicity parameter = ^aa {k, t —>■ cx)) for the active particle 
(dashed line) is slightly slower than = <hpp {k, t —)■ oo) for the passive particle, indicating 
that the glassy part formed by the active component is softer than the passive part. With 
increasing A“, the plateau heights of both correlators decrease. Above some threshold value 
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Figure 1: (Color online) (a) Dynamic scattering functions ^pp{k,t) for passive particles (solid 
lines) and <l>aa {k, t) for active particles (dashed lines) evaluated at k = km for different values of 
the activity parameter A“. The total volume fraction is (^ = 0.524 and the number fraction of active 
particles is Xa = 0.5. (b) Dependences of the non-ergodic parameters evaluated at k = km on 
the total volume fraction ip for different values of A“ and Xa = 0.5. The right inset shows for 
several values of A“ close to 1.0 for clarity. The left inset plots the dependence of ipc on A“ where 
a turnover phenomenon is apparent, wherein the dashed line refers to (p = 0.618 where the system 
shows reentrance behavior. See the text. 

of A“(about 1.044 here), the correlators hnally decay to zero and the system becomes fluid, 
wherein the active component relaxes faster than the passive one. 

Above results clearly demonstrate that increasing activity of the active component can 
shift GT to higher volume fractions. In Fig{T](b), we plot the nonergodic parameters fl 
evaluated at k = km as functions of the total fraction ip for several different values of A“. As 
it should, takes a hnite value at a certain discontinuous glass transition point (pc before 
which /I* = 0. Clearly, the value of ipc increases fast with A“ at hrst, e.g., ipc — 0.515 for 
A“ = 1.0 and (pc — 0.622 for A“ = 3.0. Strikingly, however, with further increasing of A“ to 
larger values like 5.0, (pc decreases again to a smaller value around 0.620. In the left inset of 
Fig{T](b), the dependence of ipc on A“ is shown, where the turnover phenomenon is apparent. 

The above hndings indicate an interesting type of reentrance behavior in active-passive 
mixtures for a hxed total volume fraction ip. For instance, for ip = 0.618 as indicated by 
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the dashed line in the left inset of Fig{T]^b), the system changes hrst from glass to liqnid 
and then to glass again as A“ increases. While the hrst transition from glass to liquid is as 
expected because activity can push the glass transition to higher volume fraction, the second 
one from liquid to glass again is rather counterintuitive. To get more insight, we have plot 
the nonergodic parameters and for all wave vectors in FigJU^a) for (p = 0.618 and A“ 
from 1.05 to 6.0. For A“ = 1.05 where the system is in glassy state, we see that the prohles 
of/^ and are nearly the same, with only slight differences in their values. For A“ = 2.0, 
the system is still glass, with the active part much softer than the passive part. For A“ = 3.0 
and 4.0, the system is in the liquid phase because all the /^s are now zero. Nevertheless, for 
further larger values of A“, say 5.0 and 6.0, the system becomes glassy again since both 
and have nonzero values in some k range. Interestingly, we hnd that the prohles of 
and /I* are quite different in this latter case. fl remains a relatively large value and contains 
many peaks indicating the passive component is frozen at all length scales but with a certain 
type of structure. However, /“ is only apparently nonzero for small k and is almost zero 
for k larger than some threshold value, suggesting that the active part is still liquid-like in 
short length scales but being frozen in large scales. 

We note here that a few interesting reentrance behaviors regarding glass transition have 
been reported in the literatures. For instance, for passive hard-spheres with short-range 
attractions, it was found that increasing the attraction can melt glass, however, further 
increasing attraction can lead to glass again j^. The reentrance was due to the existence 
of two qualitatively different glassy states: one with structural arrest due to caging and the 
other with arrest due to bonding. For passive hard sphere mixtures with very disparate sizes. 


reentrance behavior 


as well as multiple glasses |29| have also been reported. Reentrance 


glass transition has also been found for fluids in porous media, where for large volume fraction 
of immobile matrix particles, increasing or decreasing the number fraction of fluid particles 
may both lead to glassy states 30|, 1^. Multiple reentrant glass transitions were also found 
in a conhned system of hard spheres between two parallel walls 321] . Here we report a new 
type of reentrance glass transition behavior induced by particle activity. As demonstrated 
in FigJJl the reentrance here may be also due to existence of two types of nonequilibrium 
glassy states. For small activity, the dynamic arrest is due to the caging effect of both 
passive and active particles. Increasing activity acts as ’noise’ which can destroy the cages 
and melt the glass. If activity is too large, however, local phase separationjs-Q] of the active 
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Figure 2: (Color Online) (a)Non-ergodic parameters (solid lines) and (dashed lines) in the 
whole k range for different values of A“ at fixed total volume fraction (p = 0.618 and Xa = 0.5 
where reentrance behavior can be observed. For relatively large A“ = 5.0 and 6.0, is nearly zero 
for large k and apparently nonzero for small k. (b) Dependence of the glass transition point ipc on 
particle activity A“ for different Xa- Reentrance behavior disappears if Xa is too large. 

particles could happen such that active particles can form local clusters, which may also 
leads to dynamic arrest. In this case, the non-clustered active particles can still move like 
liquids within small length scales, such that is nearly zero for large k. It is interesting 
to check such a scenario by large scale computer simulations, which dehnitely deserves a 
separate future work. We would like to mention here that in a recent simulation work of 
active-passive mixtures j^, the authors demonstrated that the active component triggers 
phase separation into a dense and a dilute phase, where active-passive segregation were 
further observed in the dense phase with rafts of passive particles in a sea of active particles. 
We suggest that the reentrant glassy state for large activity we found in the present work 
may look like this type of active-passive segregation in dense phase. 

We now consider how the above results depend on the number fraction Xa of the active 
component. In Figj2](b), the dependences of (pc on A“ for different values of Xa are shown. 
Clearly, the reentrance behavior can occur in a large range of values of Xa- But for too 
large Xa, p>c will become saturated to a value about 0.65, which is consistent with the result 
for a pure active systemj^. Thus one conclude that the interesting reentrance behavior 
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reported above is a specific feature of active-passive mixture system and cannot occur for a 
pure active system. The data also indicates that (pc increases monotonically with Xa for a 
fixed A“, which is in consistent with the simulation work of R. Ni etal. that doping active 
particles can help crystallizing hard spheres by melting glassjw]. Note that the reentrance 
behavior does not exist with variation of Xa- 


IV. CONCLUSION 

In conclusion, we have developed a mode-coupling theory starting from Smoluchowski 
equations to study the nonequilibrium glassy dynamics in mixtures of active-passive par¬ 
ticles. This microscopic theory makes it convenient to study the relaxation of the density 
correlators and make predictions about the liquid-glass transition boundaries. In particu¬ 
lar, we have applied our theory to investigate the glass dynamics of a binary mixture of 
mono-disperse active-passive hard spheres. The theory clearly demonstrates that doping 
with active particles will push the critical volume fraction p>c for glass transition to higher 
values. In addition, we find an interesting type of reentrance behavior induced by particle 
activity. For a certain given total volume fraction where the system is in glassy state without 
activity, increasing the activity level will first melt the glass, while further increasing activity 
may lead to a glassy state again. The /c-dependent nonergodic parameters for the active 
and passive components share similar profiles in the former glass state, while they are quite 
different in the latter one, wherein the active particles are frozen in large scales but remain 
fluid-like in small ones. Such a reentrance behavior is found to be a specific feature of a 
mixture system, while it is absent in a pure active system. We believe that our present work 
can offer more perspectives in the study of collective dynamics of active-passive mixtures, 
as well as the frontier topic regarding nonequilirbium glass transition. 
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